Adiabatic pumping in a Superconductor-Normal-Superconductor weak link 
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We present a formalism to study adiabatic pumping through a superconductor - normal - su- 
perconductor weak link. At zero temperature, the pumped charge is related to the Berry phase 
accumulated, in a pumping cycle, by the Andreev bound states. We analyze in detail the case 
when the normal region is short compared to the superconducting coherence length. The pumped 
charge turns out to be an even function of the superconducting phase difference. Hence, it can be 
distinguished from the charge transferred due to the Standard Josephson effect. 

PACS numbers: 73.23.-b, 74.45. +c 
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In a mesoscopic conductor a dc charge current can be 
obtained, in the absence of applied voltages, by cycling 
in time two parameters which characterize the system 
P, 0] . This transport mechanism is called pumping. If 
the time scale over which the time-dependent parame- 
ters vary is large compared to the typical electron dwell 
time of the system, then pumping is adiabatic, and the 
pumped charge does not depcnd on the detailed timing 
of the cycle, but only on its geometrical properties. Dif- 
ferent formulations have been developed to describe adia- 
batic ramping, for example, based on scattering theory in 
Refs. HEfOl or Green's function methods in Refs. 
In the scattering approach the pumped charge per cycle 
can be expressed in terms of derivatives of the scatter- 
ing amplitudes with respect to the pumping parameters 
(Brouwer's formula Q). This result is based on the so- 
called emissivities of the system |8|, which express the 
charge that flows from a lead in response to the vari- 
ation of one parameter in the scattering region. This 
formulation requires the presence of terminals which pro- 
vide propagating channels. The scattering approach has 
been later extended to hybrid systems containing super- 
conducting (S) terminals. In Refs. 0, 0] a two termi- 
nal structure comprising one superconducting lead was 
considered. Subsequently, this approach was generalized 
to multiple-superconducting-lead systems, where at least 
one normal lead is present The presence of a nor- 

mal lead is essential for generalizing Brouwer's formula 
to these hybrid structures. 

If only superconducting leads are present, at low 
enough temperature, pumping is due to the adiabatic 
transport of Cooper pairs. Besides the dependence of 
the pumped charge on the cycle, in the superconduct- 
ing pumps there is a dependence on the superconduct- 
ing phase difference(s) (the overall process is coherent). 
Moreover, in addition to Cooper-pair pumping, in equi- 
librium, there is a contribution due to the Josephson ef- 
fect if the phase difference between the two supercon- 
ducting leads is different from zero. Up to now adiabatic 
Cooper pair pumping was studied only in the Coulomb 
Blockade regime for all-superconducting systems (super- 
conducting islands weakly connected to superconducting 



leads) [3 H H IH. In this Letter we would like to 
(partially) fiil this gap, and consider adiabatic pump- 
ing between two superconducting terminals connected 
through a normal (N) region. We focus on the regime 
of an open structure (SNS weak link), where charging 
effects are negligible. This is relevant when the normal 
region is, for example, a chaotic cavity as those used for 
normal, electronic pumps [2(. 

The derivation of a formula for the pumped charge 
makes use of the connection between Berry's phase |l6( 
and the pumped charge which we prove to be 

valid also for the SNS weak link. The resulting expression 
for the charge pumped in a period can be written in terms 
of derivatives of the Andreev-bound-statc wavefunctions 
with respect to the pumping parameters. We point out 
that there is a elose analogy of the problem studied here 
with that of pumping in a Aharonov-Bohm ring |20| . 

The system under investigation (depieted in Fig. ^) 
consists of a SNS junetion, with the weak link occupy- 
ing the region —W/2 < x < W/2. The superconducting 
order parameter is given by Aoe -1 ^/ 2 and Aoe tip ^ 2 for 
the superconductor on the left-hand-side and right-hand- 
side, respectively. The properties of the weak link can be 
varied, for example, by realizing it with a semiconduetor 
and operating on two independent external gates, indi- 
cated by Xi(t) and X 2 (t) in the figure. 

The state of such a hybrid structure \4>(t)} is the solu- 
tion of the time-dependent Bogoliubov-de Gennes equa- 
tion: 



ihd t \m) = H(t)\m), 

where the Hamiltonian 
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A(r)e^ 
V 2 - U{r, t) 



(2) 

depends on time through the two parameters: H (t) = 
H[Xi(t),X 2 (t)]. In Eq. © U(r,t) is the potential that 
takes into account the effect of the time-varying external 
gate voltages, <j){r) = cp/2[6(x - W/2) ~ Q{-x - W/2)], 
A(r) = A /2[9(x - W/2) + 0{-x - W/2)] and /i is the 
superconductor chemical potential (equal for the two S 
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FIG. 1: (Color on line.) Schematic setup for the SNS pump. 
The weak link is composed of a scattering region contacted 
to the S terminals (in orange) through two fictitious N leads 
(in light blue). The transport properties of the scattering 
region (characterized by the matrix S0) are cycled in time 
via two external gates which control the parameters X\ and 
X2- The superconducting order parameters of the S terminals 
on the left-hand-side and on the right-hand-side have equal 
magnitude Ao, but differ by a phase ip. 



leads). We now assume that the state evolves adi- 

abatically and that at any time t it is in an instantancous 
eigenstate of the Hamiltonian. The instantaneous Solu- 
tions are defined by the equation: 



H(t)\jj(t))=e(t)\i,(t)), 



(3) 



whereby t plays the role of a parameter. After a cycle of 
period r, the states returns to the initial one, but with 
an added phase factor $: 



|#r))=e^(0)>, 



(4) 



The phase $ contains both a geometrical (Berry's) and 
a dynamical contribution $ = 73 — 7d- The dynamical 
phase is simply given by 70 = 4 /J" dt (ip\H\^) 

For the SNS system, where Andreev bound states are 
formed, the condition of validity of the adiabatic approx- 
imation is that the frequency hio of the time-dependent 
parameters be much smaller than the energy difference 
between any pair of Andreev bound states, or between 
any Andreev bound states and the continuum of states 
above the gap. This implies that the pumping frequcncy 
needs at least to be smaller than the superconducting 
gap A . 

It is possible to show explicitly that the charge cur- 
rent J c h carried by a Bogoliubov-de Gennes eigenstate 
\ip) is given by the expectation value of the derivative of 
the Hamiltonian H with respect to the superconducting 
phase difference pl| : 



T 2e il\ dH \l\ 



(•5) 



The charge transferred per cycle is then defined as Q = 
f o Jch(t) dt. By assuming adiabatic evolution of the state, 
and making use of Eq. (J5J , the following relation between 



the accumulated phase and the charge transferred in a 
cycle can be written: 



d 



(0) 



The first term corresponds to the instantaneous Joseph- 
son current integrated over one period, while the second 
represents the pumped charge. Using Green's theorem, 
7b can be written in terms of derivatives with respect to 
the pumping parameter of the instantaneous eigenfunc- 
tions: 



7b = -2 / dX l dX 2 Im 
's 



dtp 
'dX~ 1 



dtp 



dX 2 



(7) 



S being the area in the parameter space spanned by 
the parameters over one cycle. In Eq. J7J the notation 
(■ |-) stands for a space integration defined by (A \B) = 
J df (f) B (r) , A and B being vectors in the Nambu 
space. 

In the short junction limit (i. e. when the distance 
W between the two superconducting interfaces is much 
smaller than the superconducting coherence length) only 
the superconducting regions contribute to the space in- 
tegration in Eq. (JJJ. The instantancous cigcnfunction, 
corresponding to the Andreev-bound-state energy ej, in 
the superconducting regions can be written as: 



with 



^s, Lj (r) = £ (».,;„ v • + h Ul y" ') Xn (y, z) 



and 



(9) 
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S, R, j 



'K 



h Rn d e 



-ik 



n ' jx ) Xn(V,z), 
(10) 

n being the transverse channel index relative to the trans- 
verse wavefunction Xn{y, z) with transverse subband en- 
ergies E n j2^ . The index j labels the different Andreev 
bound states. In Eqs. JSJ and (|TU|l k n j are particle(hole)- 
like quasiparticle wavevectors given, in the Andreev ap- 
proximation (Aq <C /i — E n for any E n ), by ■ — 



a;-: 



with k Fn - Jp(fi-E n ). The 



Nambu-space vectors b^ n ^ can be calculated with the 
following procedure: i) the eigenfunction in the fictitious 
leads in the normal regions adjacent to the supercon- 
ducting interface (see Fig. QJ are calculated along the 
lines of Ref. [23; ii) the wavefunction in the supercon- 
ductor is obtained by imposing the continuity equations 
at the interfaces within the Andreev approximation; iii) 

$S,v,j ) = 1 



the normalization condition 53i/=l r ^Vs.^.j 
is imposed. 



3 



Making use of the Andreev approximation, the pumped charge reduces to 



0p = 4e ~k L dXidX2 £ 



n, J 



1 



2Imfe 



-Im 



dX 1 dX 2 
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-Re 



*0t 



' 8X 2 dXl 



db 



/3t 



^ 

9Xi """J <9X 2 



(11) 



where the sum over j runs over the Andreev bound 
states. This is the central result of this Letter, and a 
few comments are in order. We have succeeded in ex- 
pressing the pumped charge as a function of the instan- 
taneous Andreev-bound-state eigenfunction. The vectors 
b L(R)n j depend only on the parameters of the system, 
such as the normal region scattering matrix Sq, the su- 
perconducting gap, and the superconducting phase dif- 
ference. It is clear that the pumped current can be 
written in terms of the elements of the normal-region 
scattering matrbc Sq. Equation Ijlll) is a zero temper- 
ature result, and contains only the contribution to the 
pumped charge due to the Andreev bound states. At 
finite temperatures, but still smaller than the gap, the 
contributions of the different Andreev bound states ej 
are weighted by the thermal occupation 1 — 2/(ej), being 
/ the Fermi function. At temperatures of the order of the 
gap, there is an additional contribution, not contained in 
Eq. due to the propagating quasi particles with 

energies above the gap. When superconductivity is sup- 
pressed only the latter contribution, which is described 
by Brouwer's formula, remains, leading to the usual re- 
sult for the pumped charge through a normal region con- 
nected to normal leads. 

Now let us consider the following single-channel 
parametrization for the normal-conductor scattering ma- 
tra 



Sn 



iy/9 



(12) 



choosing g and a as pumping fields (X\ = g, and X 2 — 
a). The instantaneous Andreev-bound-state energy eo i s 
simply related to the transmission probability g by 



e (t) = A 0A /l- 5 (<)sin 2 (|), 



(13) 



so that td = l/h /J" eo(t). 

The charge transfered due to the Josephson current (in 
the following named also Josephson charge) reads 



)tos = — 2e — - f dt 

h J 2 



g(t) sin(f)cos(f) 



l- 5 (i)sin 2 (f) 



(14) 



Notice that it depends on the pumping frequency ui. On 
the other hand, the pumped charge does not depend on 



frequency, but only on the geometrical prop- 
pumping cycle, and it reads 



the pumping frequency, but only on t. 
erties of the pumping cycle, and it r< 

f d 
Q p = 4e / dgda——Im 

J s d<P 



'[w 




\dg 





(15) 



Interestingly, the integrand of Eq. I|15|) turns out to be 
independent of a. 

We now consider the following sinusoidal pumping cy- 
cle: g (t) = g + Ag siri(uit) and a(t) — a + Aa sin(o;i + 0o)- 
In the weak pumping limit we assume that Ag /g <C 1 
so that the integrand of Eqs. (|14|l and (|15|l vary neg- 
ligibly during the cycle. As far as the frequency is 
concerned, the maximum value of %u) in order for the 
adiabatic hypothesis to hold is Titu < An — eo, with 

In = Aoyl — g sin 2 (< / 9/2). In order to compare Q p with 

Qjos, w e choose Ku> — 0.1(A — e ), eo being equal to the 
value of eo at <p = tt/2. Note that the adiabatic condi- 
tion breaks down for ip — or g = 0, when the Andreev 
bound state is at the gap boundary. Thus in our analysis 
we shall avoid small values of those variables. 

Figure[2]shows the pumped charge as a function of the 
superconducting phasc diffcrcncc ip for different values of 
g. Q p is a non-monotonous function of ip exhibiting a 
maximum at ip = ir. For comparison, in the inset of Fig. 
121 we plot the transferred charge due to the Josephson 
current, whose absolute value is larger, with respect to 
the pumped current, by a faetor of order hto/ Aq. How- 
ever, while Q p (<p>) is an even function of p, Qjos(v 5 ) is 
odd, so that a measure of [Q((p) + Q(—(f)]/2 will sin- 
gle out only the pumped contribution. The particular 
parity of Q p is due to the faet that a time-reversal op- 
cration implies not only the reversal of phase but also of 
the pumping trajectory in parameter space. It has also 
to be noticed (see Fig. |2J that the pumped charge is 
not quantized. In absence of Coulomb blockade, charge 
quantization oecurs only for very specific pumping eyeles 
(see, for example, Ref. |4|)- In addition the global super- 
conducting phase-cohercncc of the system, which leads 
to wave funetions extending in the two superconducting 
electrodes, also hinders charge quantization. 

To study the effect of an external magnetic field in the 
normal region, we change slightly the parametrization of 
So introducing a phase faetor exp ±i(3 in the transmis- 
sion amplitudes. As a result both the pumped and the 
Josephson charge are shifted along the <p axis by 2/3, i. e. 
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FIG. 2: Pumped charge, in units of 2e, as a function of su- 
perconducting phase difference ip computed for a sinusoidal 
weak-pumping cycle. Different lines refer to different values of 
average transmission: solid line g = 0.5; dashed line g = 0.7; 
and dot-dashed line g — 0.9. For the pumping cycl e we hav e 
chosen Ag = 0.1, Aa = 2tt and fiu = 0.1A (1 - \J\ - g/2). 
In the inset the Josephson charge is plotted as a function of 
tp for the same parameters of the main panel. 



Q p /jos{ ( p) — * Jp/josip — 2/3). For example, the maximum 
pumped charge is now reached for <p — 2(3 — -k. 

As far as detection is concerned, we wish to stress that 
for realistic parameters, using Al as superconduetor, one 
can attain sizable pumped currents of the order of 5 nA. 
A sensitive setup to currents of even smaller size can be 
realized by inserting the SNS pump in a arm of a SQUID. 

In conclusion, we have presented a formalism to study 
adiabatic charge pumping in a SNS weak link. The 
pumped charge is related to the Berry's phase accumu- 
lated in a pumping cycle by the Andreev bound-state 
wavefunctions, which can be written as a function of the 
scattering matrbc of the normal region. In the short junc- 
tion limit, the pumped charge is even with respect to the 
supercondueting phase difference. Hence, it can be easily 
distinguished from the charge transferred by the Joseph- 
son current. 
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